Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Strong lensing optical depths in a ACDM universe 



Stefan Hilbert, 1 * Simon D.M. White, 1 Jan Hartlap, 2 and Peter Schneider 2 

1 Max- Planck- Institut fur Astrophysik,Kari-Schwarzschild-Straj3e 1, D-85741, Garching, Germany 

2 Argelander-Institut fur Astronomie, Auf dem Hiigel 71, D-53121 Bonn, Germany 



o 
o 

(N 
O 

Q 

m 

(N 
> 
m 
O 
oo 
m 
o 
l> 
o 



6 



X 



5 February 2008 



ABSTRACT 

We investigate strong gravitational lensing in the concordance ACDM cosmology by 
carrying out ray-tracing along past light cones through the Millennium Simulation, 
the largest simulation of cosmic structure formation ever carried out. We extend pre- 
vious ray-tracing methods in order to take full advantage of the large volume and the 
excellent spatial and mass resolution of the simulation. As a function of source red- 
shift we evaluate the probability that an image will be highly magnified, will be highly 
elongated or will be one of a set of multiple images. We show that such strong lensing 
events can almost always be traced to a single dominant lensing object and we study 
the mass and redshift distribution of these primary lenses. We fit analytic models to 
the simulated dark halos in order to study how our optical depth measurements are 
affected by the limited resolution of the simulation and of the lensing planes that we 
construct from it. We conclude that such effects lead us to underestimate total strong- 
lensing cross sections by about 15%. This is smaller than the effects expected from our 
neglect of the baryonic components of galaxies. Finally we investigate whether strong 
lensing is enhanced by material in front of or behind the primary lens. Although strong 
lensing lines-of-sight are indeed biased towards higher than average mean densities, 
this additional matter typically contributes only a few percent of the total surface 
density. 

Key words: gravitational lensing - dark matter - large-scale structure of the Universe 
- cosmology: theory - methods: numerical 



1 INTRODUCTION 

Gravitational lensing was first discovered through strong 
lensing effects which can produc e multiple images of 
distant quasars l|Walsh et alj Il979l ) and highly distorted 
images of distant extended objects such as galax- 



ies (iLvnds fe Petrosiar] 19861: ISoucail et al.1 Il987h and ra- 
dio sources ( Hewitt et al.l 19881 ). Such effects occur when 
the surface mass density of an individual object (the 'lens') 
is comparable to that across the Universe as a whole, and 
as a result they are generated only by the most massive 
and most concentrated structures. In contrast, weak lens- 
ing, detected through the small but coherent distortion 
of the imag es of distant gala xies in the same direction 
on the sky l|Tvson et al.l ll99Cll ). is sensitive to the abun- 
dance and structure of typical nonlinear objects, so-called 
dark matter halos, and is now beginning to measure the 
statistics of th e cosmic mass distribu t ion also in the quasi- 
linear regime ([Semboloni et al.l 20061: Hoekstra et al.l 120061 : 



iHetterscheidt et al.ll2007l : iMassev et al.ll20oflT Thus gravi- 
tational lensing complements microwave background, large- 
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scale structure, and Lya forest studies, which provide infor- 
mation primarily in the quasilinear and linear r egimes (for 
recen t reviews on stron g and weak lensing, see iKochanekl 
2006: ISchneiderll2006bl) . Combining all these measures to 
constrain theories for the origin of structure requires a reli- 
able model for the nonlinear phases of evolution. 

The current standard model of cosmological struc- 
ture formation is based on cold dark matter and a cos- 
mological constant. This ACDM model has been shown 
to fit a wide variety of observations of galaxies and their 
dark halos, of galaxy clusters and galaxy clustering, of 
the structure of the intergalactic medium out to redshift 
6, and, most notably, of the detailed pattern of temper- 
ature fluctuations in the cosmic microwave background. 
The par ameters of the m odel are already highly con- 
strained (|Spergel et al.ll2007l ). Early predictions for its lens- 
ing properties wer e based on analytic models for nonlin- 
ear structure (e.g. [T urner et al. 1984; ISubramanian et al.l 
1987; Nar avan fc White] 1 19881 ). bu t reliable predictions re- 



quire numerical simulation (e.g. Bartelmann et all 1 19981 : 



I Jain et al.1 |2000| ; IWambsganss et all |200^ U~ Recent grav- 
itational lensing work has confirmed the dark halo 
structure predicted by these simulations for galaxies 
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dHoekstra et al Jl2002fc ISeliak et aLll2005l ; iMandelbaum et al 



2006 bl; ISimon et al 



200 



and clusters dComerford et al 



20061 : IMandelbaum et al.N2006al : iNataraian et all l2007h ~as 



well as for the ensemble prop e rties of the dark hal o pop- 
ulation l|Semboloni et al' 20061; iHoekstra et al.ll2006r i. Sub- 
stantial efforts are currently underway to improve all these 
measurements, and these will need to be matched by a com- 
parable improvement in the theoretical predictions. 

Additional tests and constraints may be obtained from 
observations of gravitational lensing effects. For example, 
foreground matter inhomogeneities may (de-) magnify im- 
ages of distant sources thus changing their apparent mag- 
nitude. Although expected to be small for current type 



la supernova samples (Wambseanss ct al] 19971; Holz 


1998 


Riess et al.1 1 19981; IValageas 200d: 


Amanullah et al.l 


2003 


Knop et al. 2003: Barris et alj|2004 


Riess et al. 2004). 


some 



evidence for lensing effects on the observed luminosi ty dis- 
tribution may be present in higher- redshift samples (|Wand 
2005i; |jonsson etal . 2006). For future supernova surveys, one 
should be a ble to detect such effects with higher statistical 
signif i cance dMetcalJI 19991; iMetcalf fc S ilk 1999; lMintv et al.l 
|2002| ; iMunshi fc Valageasl l2006h . This will then provide a 
further test of the standard structure formation model. 

Sufficiently massive and concentrated structures along 
the line-of-sight can give rise to multiple images, to 
strongly magnified images and to highly distorted im- 
ages, so-called giant arcs. The number of such strong- 
lensing events depends on the abundance of massive ob- 
jects and on their detailed internal structure, both of which 
are sensitive to the background cosmology. Currently a 
much-debated question is whe t her t h e observed frequency 
of gi a nt arcs |Luppino et alj 1 19991 ; IZaritskv fc Gonzalez! 
2003; iGladders et all I2003D is too high compared to 
predictions based on ACDM models with parameters 
favoured by other observations (iBartelmann et al.l [l998 ; 



Meneghetti et al.ll2000l; lOguri et al.ll2003l; iDalal et all 2004 



Wambseanss et alj|2004l; iLi et alT 20051; iHoresh et al 



2005 



Li et alj|2006l : IWu fc Chiuehll2006l ; iMeneghetti et al.ll2007h 
The problem appears particularly pressing at higher red- 
shift, but available simulation results are not good enough 
to establish a clear discrepancy. 

In this paper, we will study the optical depth for a va- 
riety of strong lensing effects as a function of source redshift 
in the standard ACDM cosmology. In particular, we will es- 
timate the fraction of images that are highly magnified by 
gravitational lensing, that have a large length-to-width ra- 
tio, or that belong to multiply imaged sources. In addition, 
we compare the effect of foreground and background mat- 
ter to that of the primary lens in generating these optical 
depths. 

The results presented here were obtained by shooting 
random rays through a ser ies of lens planes cr eated from 
the Millennium Simulation l|Springel et alj|2005h . This very 
large A-body simulation of cosmic structure formation cov- 
ers a volume comparable to the largest current surveys 
with substantially better resolution than previous simula- 
tions used for ray-tracing studies. Our set of lens planes 
represents the entire mass distribution between source and 
observer, allowing us to quantify the effects of foreground 
and background matter. In addition we have ensured that 
our ray-tracing techniques take full advantage both of the 
statistical power offered by the large volume of the simula- 



tion, and of its high spatial and mass resolution. This allows 
us to make more precise statements about model expecta- 
tions than has previously been possible. 

Our paper is organised as follows. In Sec. [2] we de- 
scribe how we trace representative rays through the Millen- 
nium Simulation. Results for the magnification distribution 
as a function of source redshift are presented in Sec. 13.11 
Strong-lensing optical depths are then discussed in Sec. 13.21 
In Sec. 13.31 the mass and redshift distribution of the objects 
which cause strong lensing are examined, and we demon- 
strate that the errors induced by the finite volume and reso- 
lution of the simulation are relatively small. Biases induced 
by additional structure in front of or behind the principal 
lens are examined in Sec. l3.4l and are also found to be small. 
The paper concludes with a summary and outlook in Sec. U 



2 METHODS 

2.1 The multiple-lens-plane approximation 

In the multiple-lens-plane approximatio n (see, e.g., 

iBlandford fc Naravanl [l986l ; ISchneider et"ai1ll992T ), a finite 
number of planes are introduced along the line of sight, 
onto which the matter inhomogeneities in the backward 
light cone of the observer are projected. Between these 
lens planes, light is assumed to travel on straight lines. 
Light rays are deflected only when passing through a lens 
plane. The deflection angles may be calculated from the 
gradient of a lensing potential, which is connected to the 
projected matter distribution on the lens planes via a 
Poisson equation. 

A light ray reaching the observer from a given angular 
position 6 can then be traced back to the angular position (3 
of its source at given redshift z s , thereby defining the lensing 
map 



1 — ► f3 



(1) 



from the image plane P 1 to the source plane P s . The distor- 
i.e. the Jacobian of the map, quantifies 



tion matrix A 



90 

de > 



the magnification and distortion of images of small sources 
induced by gravitational lensing. The (signed) magnification 
fi of an image is given by the inverse determinant of the dis- 
tortion matrix: 



H = (dot A) 1 . 
The decomposition l|Schneider et al.lll992l ) 



COS If 

- sin <p 



sm if 
cos ip 



1 — k — 71 

-72 



1 



-72 

K + 71 



(2) 



(3) 



of the distortion matrix defines the image rotation angle tp, 
the convergence k, and the complex shear 7 = 71 + i72- The 
reduced shear g = 7/(1 — k) determines the major-to-minor 
axis ratio 



i + lsl 



(4) 



of the elliptical images of sufficiently small circular sources. 
The determinant and trac e of the distortion mat rix may be 
used to categorise images l|Schneider et al.iri992l ')FI 

1 In lSchneider et al.l lll992t) . the trace and determinant of the 
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• type I: det A > and tr A > 0, 

• type II: det A < 0, 

• type III: det A > and tr A < 0. 

In all situations relevant for this work, images of type II 
and type III belong to sources that have multiple images. In 
the following, we will often consider type II and III images 
together as 



• type II V III: det A < or tr A < 0. 



In this paper, we want to study how often one can ex- 
pect to observe images with certain lensing properties, e.g. 
highly magnified or strongly distorted images. In order to 
quantify the frequency of rays with a given property p, we 
define the optical depth 



where 



1p(0) = 







, Jg A 2 ei p {e) 

if ray(0) has property p, 
otherwise. 



(5) 



For a uniform distribution of images in the image plane, Tp 
estimates the fraction of images of sufficiently small sources 
that have the property p. Furthermore, we define the optical 
depth 



s 

T„ = 



jg d 2 Hfl)r%(fl) 

Ld*o He)]- 1 



(6) 



to estimate the average fraction of images with properties 
p for a uniform distribution of sources in the source plane. 
These optical depths (assumed to be smooth) will be used 
to define corresponding probability density functions (pdf) 
for the magnification: 



r i/s 



(7) 



Compared to t p , the optical depth Tp takes into ac- 
count that areas in the image plane with higher magnifica- 
tion map to smaller areas in the source plane. This aspect 
of magnification bias leads to a lower image density in ar- 
eas of higher magnification for volume-limited surveys. In 
magnitude-limited surveys, however, magnification can push 
images that would otherwise be too faint to be observed 
above the detection threshold. This aspect of magnification 
bias counteracts the previous one, but will not be discussed 
in this paper since it depends sensitively on the luminosity 
distribution of the source population. 

Note that Tp and Tp differ from the optical depths 



f pS d 2 /31 P (/3) 



(8) 



discussed, e.g., by ISchneider et al.l (|l992h , which quantify 
the fraction of sources whose images have certain proper- 
ties. The methods used in this paper do not yield enough 



symmetric part [i.e. second factor on the r.h.s. of Eq. ((3j ] of the 
distortion matrix A have been used. However, the determinant 
and the sign of the trace of A and its symmetric part are identical 
for \(p\ < it/2. 



information to deduce fp in general. We will therefore re- 
strict our discussion to Tp and Tp. 

Roughly speaking, Tp weights lensing events by their 
area on the sky, Tp weights by the number of images, and 
fp weights by the number of sources. In the absence of multi- 
ple images, Tp and would be identical. For strongly lensed 
properties such as considered in this paper, they can be sig- 
nificantly different: For Tp , each of the multiple images of 
a source contributes individually, whereas all images of the 
same source contribute a single event event to fp. Conse- 
quently, for a given number N s of uniformly distributed 
sources in the source plane, N S Tp is the expected number 
of images with property p, whereas N s fp gives the expected 
number of sources. Obviously, the number of images is easier 
to count in observations than the number of sources. 



2.2 The lens planes 

Our methods for reconstructing the observer's backward 
light cone, for splitting it into a series of lens planes, and 
for calculating the matter distribution and deflection an- 
gles on these pla nes are generally similar to those used by 
Ijain et all |2000t) . They differ, however, in a number of im- 
portant details which reflect our wish to take full advantage 
of the unprecedented statistical power offered by the large 
volume and high spatial and mass resolution of the Millen- 
nium Simulation. Here, we give a brief outli ne of our al- 
gorith ms, reserving a detailed description for iHilbert et al.l 
(|2007h . 

The Millennium Simulation i|Springel et al.ll2005l ) is an 
iV-body simulation of cosmological structure formation in a 
flat ACDM universe with a matter density of Q,m = 0.25 
(in terms of the critical density), a cosmological constant 
with Qa = 0.75, a Hubble constant h — 0.73 in units of 
100 kms _1 Mpc _1 , a primordial spectral index n = 1 and 
a normalisation parameter erg = 0.9 for the linear density 
power spectrum. The simulation followed N « 10 10 particles 
of mass m p = 8.6 x 10 8 /i _1 Mq in a cubic region of comoving 
side L = 500/i _1 Mpc (assuming periodic boundary condi- 
tions) from redshift z = 127 to the pr esent using a TreePM 
version of GADGET-2 (|Springell |2005| 'I. The force softening 
length was chosen to be 5ft" 1 kpc comoving. Snapshots of 
the simulation were stored on disk at 64 output times spaced 
approximately logarithmically in expansion factor for z > 1 
and at roughly 200 Myr intervals after z — 1. 

Since the fundamental cube of the simulation is too 
small to trace rays back to high redshift in a single repli- 
cation, we have to make use of the periodicity to construct 
our light cones. In order to reduce the repetition of structure 
along long lines-of-sight (LOS) through this lattice-periodic 
matter distribution, we chose the LOS to be in the direc- 
tion n = (1,3,10). This results in a comoving period of 
5.24h~ Gpc along the LOS, giving the first image of the 
origin at z — 3.87. It also allows us to maintain periodicity 
perpendicular to the LOS with a rectangular unit cell of di- 
mension 1.58/i _1 Gpc x 1.66/i _1 Gpc comoving. This period- 
icity allo ws us to use Fast-Fourier- Transform (FFT) me th- 
ods (e.g.. ICoolev fc Tukevlll965l : iFrigo fc Johnson! l2005h to 
obtain the lensing potential and its derivatives on the lens 
planes. 

To construct the matter distribution in the observer's 
backward light cone, we partitioned space into a series of 
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redshift slices, each perpendicular to our chosen LOS and 
containing the part of the light cone closer to one of the 
snapshot redshifts than to its neighbours. The matter dis- 
tribution within each such slice was then approximated by 
the stored particle data at the time of the corresponding 
snapshot, was projected onto the lens plane, and was placed 
at the comoving distance corresponding to the snapshot's 
redshift. In order to reduce the shot noise from individual 
particles, we employed an adaptive smoothing scheme. Each 
particle was smeared out into a cloud with projected surface 
mass density 



f 3m p A Ix-gpl' N 2 



\X Xp I <C Tp , 



(9) 



where x denotes comoving position on the lens plane, x p is 
the projected comoving particle position, and r p denotes the 
comoving distance to the 64th nearest neighbour particle in 
three dimensions (i.e. before projection). 

To calculate the lensing potential from the projected 
matter density, we used a particle-mesh particle-mesh 
(PMPM) method. The effective spatial resolution of the Mil- 
lennium Simulation is about 5ft" 1 kpc comoving in dense re- 
gions, where the particles' smoothing lengths become com- 
parable to softening length of the simulation, i.e. r p ~ 
5ft" 1 kpc. Hence, a mesh spacing of 2.5ft" 1 kpc comoving is 
required to exploit the numerical data without degradation. 
A single mesh of this spacing covering the whole periodic 
area of the lens plane (i.e. 1.58ft" 1 Gpc x 1.66ft" 1 Gpc co- 
moving) would, however, be too demanding both to compute 
and to store. We therefore split the lensing potential <!> into 
long-range and short-range parts defined in Fourier space 
by: 



*ion g (fc) = *(fc) exp (-Z S pii t fc ) , and 
*sho rt (fc) = tt(fc) [1 - exp (-4 lit fc 2 )] 



(10a) 
(10b) 



The comoving splitting length l Bp ut = 0.175ft" 1 Mpc charac- 
terises the spatial scale of the split. We use a 16384 x 16384 
mesh covering the whole periodic area of the lens plane to 
calculate ^long from the projecte d surface mass density by 
an FFT (|Frigo fc Johnson! |20o3 ) method. The long-range 
potential is calculated once and then stored on disk for each 
of the lens planes. To calculate *& s hort, a fine mesh with 
2.5ft" 1 kpc spacing is used. This mesh only need cover a rel- 
atively small area around regions where the potential is re- 
quired, i.e. where light rays intersect the lens plane. Because 
of the short range of "fshort, periodic boundary conditions 
can be used on the fine mesh, provided points close to its 
boundary are excluded from subsequent analysis. Therefore, 
FFT methods can used without 'zero padding' to calculate 
^short on the fine mesh. The long- and short-range contri- 
butions to the deflection angles and shear matrices (i.e. the 
second derivatives of the lensing potential) are calculated on 
the two meshes by finite differencing of the potentials. The 
values between mesh points are obtained by bilinear inter- 
polation. The resulting deflection angles and shear matrices 
at ray positions can then be used to advance the rays and 
their associated distortion matrices from one plane to the 
next. 



2.3 Ray sampling 

In order to estimate optical depths and magnification dis- 
tributions, we shoot random rays through our set of lens 
planes. In doing this, we neglect correlations between the 
matter distributions on different lens planes. This allows us 
to pick random points on each lens plane as we propagate 
rays back in time, significantly simplifying and accelerating 
our code. Since the comoving separation between the lens 
planes is large (~ 100ft -1 Mpc), and the shear matrices are 
dominated by small-scale structure (< lft" 1 Mpc), this as- 
sumption is very well justified for the purposes of the current 
paper. 

On each lens plane, 40 fields of about 40ft" 1 Mpc x 
40ft" 1 Mpc comoving were selected at random. Within each 
of these fields, the shear matrix was calculated at 16 million 
random positions by our PMPM algorithm. The resulting 

6.4 x 10 8 shear matrices for each plane were then stored 
on disk together with the positions at which each had been 
computed. 

The shear matrices from all lens planes were next com- 
bined at random to produce 640 million simulated LOS from 
the observer back to high redshift. Imagining a small circular 
source at the position where each of these LOS intersects a 
particular lens plane, we can combine the shear matrices of 
the ray from all lower redshift planes to obtain the trace tr A 
of the distortion matrix, the magnification fi = (detA)" 1 , 
and the length-to-width ratio r of the source's image. The 
measured fractions of rays with certain properties, e.g. a 
large magnification, can then be used to estimate the corre- 
sponding optical depth to the redshift of the chosen plane. 

This sampling method assumes that the rays are uni- 
formly distributed in the image plane. The observed number 
N p of rays with a particular property p, when compared to 
the total number of rays TV, is thus a straightforward Monte- 
Carlo estimate (without importance sampling) of the optical 
depth Tp-. 



Np 
TV 



(11a) 



The calculation of the corresponding optical depth Tp re- 
quires using the individual magnifications fj,(i) of the rays 
i = 1, . . . , TV as statistical weights: 



Ef =1 Ia^COI 



(lib) 



The matter distribution on our lens planes is guaran- 
teed to be periodic, smooth and non-singular as a result of 
the adaptive smoothing we use. Furthermore, a large and 
random area of the image plane is mapped onto an equally 
large area in the source plane by our lensing map Q. Thus 
a representative ray sample should satisfy [see Eq. (|A2|) in 
Appendix IX] : 



(12) 



For our ray sample, we find this relation to be satisfied to 
quite high precision; for all source planes the deviation is 
smaller than 0.003. 

Our ray sampling technique neglects correlations be- 
tween the structure on different lens planes. The effects of 
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Figure 1. The probability density pdf (fi) of the magnification fi. 
(a) Logarithmic plot comparing pdf : (/x) (solid line) and pdf s (/^) 
(dashed lines) for a source redshift of z s = 2.1. (b) Linear plot 
comparing the pdf ! (/x) around fi = 1 for different source redshifts, 
2 s = 1.1 (solid line), z s = 2.1 (dashed lines) and z s = 5.7 (dotted 
lines). 

the lens environment on scales smaller than 100/i -1 Mpc co- 
moving should be correctly represented, however, and the 
effects of uncorrelated fluctuations in the density of fore- 
ground and background matter are also included properly. 
This simple procedure should thus give accurate results in 
the context where we use it, but we note that it does not al- 
low the construction of extended images of extended sources. 
This can be done by relatively straightforward extensions of 
our methods which we reserve for future papers. 



3 RESULTS 

3.1 The magnification distribution 

From the magnifications of our random rays, we have 
estimated the probability density functions pdf ! (/i) and 
pdf s (/i). These are compared for a source redshift z s =2.1 
in Fig. [T^. One can readily see the stronger fall-off to high 
magnification for pdf s (/i), which is a consequence of the 
relation pdf s (/i) « \fj,\~ l pdf(jj.) [see Eq. (|A3|l in the Ap- 
pendix] . The asymptotic behaviou r predicted from catastro- 
phe theory (|Schneider et al.lll992l ). i.e. pdf : (/i) oc /i~ 2 and 
pdf s (/i) oc fj,~ 3 , is reached for magnifications fj, > 20. 

Probability density functions pdf s (/i) for different 
source redshifts z s are shown in Fig.[lp. With increasing z s , 
the peak of the distribution broadens and moves to lower fi, 




. 9 
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(b) width of: 




1 2 3 4 5 6 

z s 

Figure 2. The peak position (a) and the width (b) of the prob- 
ability densities pdf : (^t) (solid lines) and pdf s (/x) (dashed lines) 
for the magnification fi as a function of source redshift z s . The 
distributions get broader and more skew with increasing z s . 

whereas the high-^i tail increases in amplitude. The peak po- 
sitions and the widths of the pdfs, i.e. their modes /ipeak and 
full-widths-at-half-maximum FWHM, are plotted as func- 
tions of source redshift in Fig. [2] The shift of the peak 
with increasing redshift balances the heavier tail so that 
/Atpdf S (M)a> w 1 for all redshifts [see Eq. flU)]. These 
results for pe ak position and widt h are in goo d agre ement 
with those of lValageasl (|2000h and lFluke et al.1 (|2002l ). who 
considered similar cosmologies. 

There is a lower bound to the magnification of images 
of type I, which is realized for rays which propagate through 
empty cones (i.e., for which the matter density vanishes 
along their path) and which ar e subject to no shear e f- 
fects (e.g., iDver fc Roederl 1 19721 ; ISeitz fc Schneiderl I1992T I. 
No light ray can diverge more strongly than such an empty- 
beam ray. Lower magnifications can only be produced for 
overfocussed rays which then belong to Type II or III. A 
simple way of calculating this lower bound fimin for a flat 
universe is given in Appendix [B] The steep rise in the prob- 
ability density of the magnification at /i ps 0.83 seen in 
Fig. Ola) is substantially larger than the theoretical bound 
fJ-mia = 0.69 for z s — 2.1, indicating that there are no real 
empty cones in a realistic universe, w hich is in agreement 
with findings by I Vale fc White! l|2003l ). 

3.2 Strong- lensing optical depths 

Sufficiently concentrated matter clumps between distant 
light sources and the observer can give rise to highly mag- 
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Figure 3. Optical depths for images of small circular sources 
of type II (solid lines), of type III (dash-dotted lines), with large 
magnification (dashed lines) and with large length-to- width ratio 
(dotted lines), assuming a uniform distribution of images in the 
image plane (a), and a uniform distribution of sources in the 
source plane (b). Note that the optical depths are significantly 
smaller in the latter case, and that the relative optical depths for 
different types of strong lensing are not the same in the two cases. 



nified, strongly distorted, or multiple images. We refer to 
such phenomena as strong lensing. In order to quantify the 
amount of strong lensing expected in a ACDM universe with 
the parameters of the Millennium Simulation, we used our 
large set of random rays to estimate: 

• the fraction with det A < (type II), 

• the fraction with det A > and trA < (type III), 

• the fraction with det A < or trA < 0, i.e. the sum of 
the two previous classes (type II V III), 

• the fraction with a length-to-width ratio r > 10 for im- 
ages of sufficiently small circular sources, and 

• the fraction with magnification > 10. 

The corresponding optical depths Tp(z s ) and Tp(z s ) are plot- 
ted in Fig. [3] as functions of the source redshift z s . Since all 
the image properties we consider are (either by definition or 
at least statistically) associated with large magnifications, 
the optical depths Tp (z s ) are always a factor 5 to 20 smaller 
than the corresponding Tp(z s ). 

The optical depths for r > 10 are 2 to 20 times smaller 
than those f or |/x| > 10. Si milar resu l ts hav e been found by 
iDalal et all (|2004 ) and by iLi et al.1 |2005l ). Evidently, the 
optical depth for highly magnified images does not provide 
a reliable estimate for the probability of images with a large 
length-to- width ratio. 



The optical depth tv >10 may be a reasonable approx- 
imation to the optical depth for giant arcs with length-to- 
width ratio r > 10, since both finite source size and finite 
source ellipticity affect this particular property only weakly. 
Moreover the two effect s work in opposite directions. For 
example, ILi et al.l ((2005) found that the optical depth for 
r > 10 is almost identical to the optical depth for arc im- 
ages with length-to-width-ratios > 10 of elliptical sources 
with an effective diameter of 1 arcsec. 

All our optical depths show a strong dependence 
on source redshift , simi lar t o that p r evious ly noted by 
IWambsganss et all l|2004 ) and ILi et all (|2005h . Comparing 
to their results in detail, however, we find a somewhat 
stronger redshift dependence, resulting in 2 to 3 times higher 
optical depths at z s > 0.5. There are various possible ex- 
planations for this discrepancy. The simulations they use 
have 3 to 20 times worse ma ss resolution than the Mil len- 
nium Simulation. In addition. Wambsganss et al.1 (|2004h re- 
duced their spatial resolution on lens planes with increasing 
redshift, thereby potentially missing some low mass lenses, 
whereas we always use the full spatial resolution permit- 
ted by the fo rce resolution of the Mil lennium Simulation. 
Furthermore, IWambsganss et all (|2004l ) measured the frac- 
tion of multiply imaged sources for which at least one image 
has \fi\ > 10. Thus they do not include sources with a sin- 
gle highly magnified image, as can occur for a marginally 
subcritical lens. Moreover, they do not account for the fact 
that multiply imaged sources c an give ri s e to m ore than one 
image with \fi\ > 10. Similarly. ILi et al.1 (|2005l ) only consid- 
ered the magnification /j, and length-to-width-ratio r of the 
brightest image of multiply imaged sources for the calcula- 
tion of the optical depths for > 10 and r > 10, although 
they took into account all giant-arc images for a given pop- 
ulation of extended sources in their calculation of the opti- 
cal depth for giant arcs. Furthermore, they only considered 
massive clusters and neglected possible contributions from 
foreground and background matter. These decisions caused 
them to underestimate the total cross sections for strong 
lensing (at least for smaller sources). 

Fold singularities generically produce strongly magni- 
fied and distorted images as pairs, and sources near cusps 
may giv e rise to one or three s trongly magnified or distorted 
images l|Schneider et al.lll992l ). Therefore, sources with mul- 
tiple highly magnified or strongly distorted images might 
explain a substanti al part of the discr e pancy bet ween our 
result s and those of IWambsg anss et al.l |2004) and lLi et al.l 
(2005). With our methods, we cannot quantify this effect, 
but we discuss the other effects further in the following sec- 
tions. 



3.3 Lens properties 

In most cases, the properties of strongly lensed rays, i.e. rays 
with det A < 0, tr A < 0, > 10, or r > 10, are predom- 
inantly caused by a single matter clump along the line of 
sight. We refer to this as the lens of the ray. In order to find 
these clumps and to study their properties, we determined 
for each strongly lensed ray those lens planes which were suf- 
ficient to produce the relevant property in the single-plane 
approximation. Only 2 x 10 -4 of all rays had more than one 
'sufficient' plane, and we will simply ignore these rays in the 
following. On the other hand, there was no sufficient plane 
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for up to 41 percent of the rays, depending on source redshift 
and the property considered, so we will discuss these cases in 
more detail in Sec 13.41 For most strongly lensed rays, how- 
ever, this simple criterion identifies exactly one lens plane. 
The redshift of this plane was then taken as the lens red- 
shift z L for the ray. For rays of type II V III (i.e. rays with 
detA < or trA < 0), the resulting lens redshift distribu- 
tion is illustrated in Fig. [4] Here we plot the cross section 
dr l /dz L as a function of lens redshift z L for various source 
redshifts z s . Even for sources at redshift z s — 5.7, most of 
the lenses have redshift z L < 2.5. The relatively low cross 
section at higher z L reflects both the lower abundance of 
massive halos and the less favourable geometry for lensing 
at these redshifts (see Fig. [6}. The lens redshift distribu- 
tions for rays with > 10 and with r > 10 are almost 
indistinguishable from that of Fig. [4] despite the different 
total optical depths. 

We studied not only the redshift of the clumps acting 
as strong lenses, but also their masses. All significant mat- 
ter concentrations have already been identified as DM halos 
in the simulation and their masses and central positions are 
available in the simulation archive. We first projected the 
centres of all halos onto the lens planes in the same way as 
was done for the particles. For each strongly lensed ray and 
for all DM halos close to the point where the ray intersects 
a lens plane, we determined the ratio M L /b, where M L is 
the conventional halo mass (defined as the mass within a 
sphere with mean enclosed density 200 times the cosmolog- 
ical mean), and b denotes the impact parameter of the ray 
with respect to halo centre. 

The DM halo with the largest M L /b on the sufficient 
plane was then defined to be the lens of ray. We discarded 
from further analysis those 3 to 6 percent of rays for which 
the largest value of M L jb was not at least ten times the sec- 
ond largest value. This cut removed all cases where one could 
not easily separate the influence of several neighbouring DM 
halos, for example in merging clusters[f| The resulting distri- 
butions of lens masses for rays of type II V III, with > 10, 
and with r > 10 are compared in Fig. [5^, where the cross 
sections dr l /<91ogM L are plotted for z s =2.1 as a function 
of lens halo mass M L . Although the corresponding total op- 
tical depths are quite different, their lens mass distributions 
are very similar. There is only a small shift toward lower 
masses for rays with \fi\ > 10, and r > 10 compared to rays 
of type II V III. In the following we will restrict discussion 
to the latter for simplicity. 

In Fig.(5]D, the cross section dr l /d\og A1 L is plotted as a 
function of mass for type-II V III rays and for various source 
redshifts z s . The measured cross section dr 1 /<91ogM L 
vanishes for masses below 4x 1O 12 /i _1 M0 and above 
4 x 10 15 /i _1 Mq. The main contribution to the optical depth 
r 1 comes from halos with 10 13 /i _1 M < M L < 10 15 /i _1 M . 
For higher source redshifts, the cross section has more weight 
at lower masses. 

The upper mass limit for the cross section simply 
reflects the fact that there are no halos more massive 
than 4 x 10 h~ Mq in the simulation. However, the cross 



2 The choice M L /6 is somewhat arbitrary. We also tried Af L /6 2 , 
but a different halo was chosen only in a few cases, all of which 
were removed by our imposed ratio cut. 



section decreases rapidly with increasing M L already for 
il/ L > 10 15 /i -1 Mq. The exponential decrease in halo abun- 
dance with increasing mass apparently dominates over the 
increasing cross section of individual halos. 

At all source redshifts there is a significant contribution 
from halos with M L < 10 14 /i _1 M©. For the small sources 
considered here, the set of DM halos causin g strong l ensing 
exte nds to substantiall y lower masses than iLi et all (|2005f ) 
and iDalal et all £2004) suggest for halos producing giant 
arcs. The Millennium Simulation has 12 to 20 times bet- 
ter mass resolution than the simulations used by these au- 
thors. This provides a considerably better representation of 
the central halo regions which produce strong lensing. In ad- 
dition, halos are inefficient in generating strongly distorted 
images for sources with angular extent comparable to or 
larger than their Einstein radii. Thus, neglecting or incor- 
rectly treating halos below a given mass (e.g. because of 
limited mass resolution) has a larger effect on the optical 
depths for small sources than on those for extended sources 
- especially for high source redshifts. Together these effects 
may explain why we find 2 to 3 times la rger optical dep ths 
at high redshift than the values given bv lLi et al.l (I2005T ). 

What sets the lower mass limit for the cross sec- 
tion d^/dlogM 1 ? The nominal mass resolution for 
identifying DM halos in the Millennium Simulation is 
about 10 10 /i -1 Mo, so there are plenty of halos with 
M L < 4 x 10 12 /i _1 M Q . The identification limit cannot, 
in itself, explain the lack of halos less massive than 
4 x 10 12 /i _1 M o in our sample. However, the regions capa- 
ble of causing strong lensing are very small for low-mass 
halos, so our cross section estimates may be limited by the 
resolution of our lens planes; critical regions with a diam- 
eter below the mesh spacing or the effective gravitational 
smoothing scale are not resolved. In order to estimate the 
mass limit i nduced by these effe cts, we considered spherical 
NFW halos l|Navarro et aljl997l ) with concentration param- 
eter 

, m -l L\ 9.59 ( M L V 0102 

c(M '* )= i^{ w^mJ (13) 

determined by halo mass and redshift (|Dolag et alj|2003 ). 
For given lens and source redshifts, there is a minimum 
lens mass for which the Einstein radius exceeds the reso- 
lution limit of our lens plane. We take the latter to be 5 kpc 
comoving, thus requiring a minimum of four mesh points 
across the Einstein diameter. This limit takes into account 
not just the limit imposed by the mesh spacing of 2.5h~ kpc 
(~ 1.7/i _1 kpc in radius), but also the force softening and the 
smoothing for particles in the halo cores. The solid line in 
Fig. [6] shows the resulting minimum mass M L as a function 
of z L for z s — 5.7. The shading in this plot gives the cross 
section d 2 T I /dz L dlogM Ij . The region of the (z L , log M L )- 
plane with non-zero cross section is bounded above by the 
largest halo mass at each redshift (the dashed line). The 
lower boundary of this region lies slightly below our ana- 
lytic estimate of the resolution limit (the solid line). About 
6 percent of the cross section is below the analytic estimate. 
Some deviation is expected because this estimate does not 
include the effects of intrinsic ellipticity, asymmetries and 
substructure, and of the scatter in concentration of halos of 
given mass. It also neglects scatter due to additional mat- 
ter inhomogeneities along the LOS. These effects should re- 
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suit in strong lensing by somewhat lower mass halos than 
our s i mple spherical model would indicate (pVIeneghetti et al.1 
120071 : iHennawi et~aDl2007l : iFedeli et al.ll2007l ). In particular, 
a high-concentration, prolate halo with its major axis along 
the LOS has a greatly enhanced cross section relative to a 
spherical halo of the same mass with a typical concentra- 
tion. Moreover, ellipticity and scatter in concentration lead 
to larger cross sections on average compared to spherical 
NFW halos of the mean concentration. 

Besides the resolution limit due to the mesh spacing 
and local smoothing scale on the lens plane, there is another 
factor limiting the resolution: The critical regions of the rel- 
evant halos may not have been simulated to the accuracy 
required to get fully converged results in the face of dis- 
creteness effects resulting from the relatively small number 
of particles in t h ese re gions. According to the criteria given 
by iPower et al.l l|2003l ). only the most massive halos in the 
Millennium Simulation at 0.3 < z < 2 should have spheri- 
cally averaged density profiles which are fully converged at 
radii comparable to their Einstein radius. At lower masses 
and at other redshifts the particle nu mber in the inner re - 
gions is below the value advocated by IPower et al.l |2003). 
(In contrast, the softening length of the Millennium Sim- 
ulation appears adequate to avoi d major prob l ems.) When 
the particle number is too small, IPower et al.l l|2003h show 
that simulations typically underestimate the central concen- 
tration of a halo, implying a reduction in its strong lens- 
ing cross section. This effect is a relatively slow function of 
particle mass. In addition, strong lensing depends on the 
projected de nsity distribution r ather than the 3-D density 
profile which IPower et all (|2003l ) studied; there are typically 
two to three orders of magnitude more particles projected 
within a halo's Einstein radius than there are within a cen- 
tral sphere of radius r_g. Thus it is unclear how seriously the 
under-resolution of halo cores will affect the cross sections 
we calculate. 

To obtain a rough estimate of how much optical depth 
we lose due to resolution effects, we calculated the cross 
section chn vm /<9 log M L approximating all DM halos in the 
Millennium Simulation by spherical NFW halos while either 
(i) taking into account or (ii) disregarding halos with an Ein- 
stein radius te < 5/t -1 kpc comoving. In doing this, we used 
the measured maximal circular velocity of each halo to esti- 
mate its concentration parameter, rather than assuming the 
concentration to be given by Eq. 1)131). The cross sections 
obtained for the two cases are shown in Fig. [7] Due to the 
scatter in the halo concentrations, the analytical estimate 
excluding halos with Einstein radius tb < 5/t _1 kpc extends 
to lower masses than the limit calculated above by assum- 
ing Eq. (|13|l for all halos. The estimate for Tn VII i including 
halos with Einstein radius te < 5/i _1 kpc is only about 15 
percent larger than the estimate excluding such halos. More- 
over, there is no significant contribution to the full estimate 
from halos below a few times 1O 12 /i -1 M0. Indeed, a de- 
tailed analysis shows that the strong-lensing cross section 
of spherical NFW halos with M L < W^h' 1 M decreases 
exponentially with decreasing mass, and is not compensated 
by the increasing number of halos. 

Fig. [7]also shows the two spherical halo-based estimates 
scaled up by a factor of 3. The curve neglecting halos with 
small Einstein radii is then a good match to the histogram 
derived directly from the simulation. Hence, the cross sec- 
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Figure 4. The cross section dr l /dz L for rays of type II V III as a 
function of lens redshift z L for sources at redshift z s = 1.1 (solid 
line), z s = 2.1 (dashed line), and z s = 5.7 (dotted line). Even for 
high redshift sources, the typical lens redshift is relatively low. 



tion cVnvnj/dlog M L is about three times as large for sim- 
ulated halos as for spherical NFW halos, at least for halo 
masses M L > 4 x W 13 ^ 1 M where the resolution limit 
is unimportant. If this result applies also at lower mass, it 
implies that limited resolution does not lead us to under- 
estimate total optical depths substantially, perhaps only by 
about 15%. Note that the missing cross section corresponds 
to very small image splittings A8 1 < 1 arcsec, a scale where 
the gravitational effects of the baryons in the central galaxy 
are expected to be important. 



3.4 Effects of additional matter along the line of 
sight 

For some rays with det A < 0, tr A < 0, \fi\ > 10, or r > 10, 
there is no individual lens plane that is sufficient to produce 
the relevant property in the single-plane approximation. The 
fraction of rays for which this is the case is shown in Fig. [8] 
It increases with increasing source redshift, is lowest for rays 
of type IIVlII, and is highest for rays with r > 10. This frac- 
tion gives an indication of the extent to which foreground 
and background material affects the strong lensing optical 
depths we have estimated. Such material is expected to have 
no or little effect on average for the image(s) of a randomly 
chosen object of given redshift. However, by selecting rays in 
the extreme tail of lensing distributions, we may be signifi- 
cantly biased towards lines-of-sight for which the additional 
material enhances the effect of the primary lens. Fig. [8] sug- 
gests that additional material is particularly effective in en- 
hancing the probability of highly distorted images (e.g. giant 
arcs), presumably because these are sensitive to the lensing 
map in a narrow region around its critical lines. One should, 
however, bear in mind that not only do certain directions 
gain a considered property through the primary lens being 
supplemented by additional LOS material, but other regions 
lose the same property because the primary lens is counter- 
acted by lower than average additional material. Therefore, 
the fractions of Fig. [8] do not reflect the overall contribution 
of line-of-sight material to our cross sections. 

In general, the effects of foreground and background 
material are relatively weak. In the cases where there is 
no single plane which is sufficient to generate the relevant 
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Figure 5. The cross section 9r I /91ogAf L as a function of the 
mass of the lensing halo M L (see text). Panel (a) compares rays 
of type II V III (solid line), rays with > 10 (dashed line), and 
rays with |r| > 10 (dotted line) for source redshift z s = 2.1. (b) 
compares rays of type II V III for sources at redshift z s = 1.1 
(solid line), z s = 2.1 (dashed line), and z s = 5.7 (dotted line). 
The lens mass distribution is almost independent of the type of 
strong lensing event, but it shifts towards lower masses for higher 
source redshift. 



property, there is still usually a single plane dominating the 
lensing effects. As an example, we determined for each ray 
with \fi\ > 10 and z s = 5.7 the planes which gave rise 
to the largest and second largest magnifications, ^i st and 
H2nd, resp., in the single-plane approximation. The cumula- 
tive distribution of these magnifications is shown in Fig. [5] 
Even though the fraction of rays with /ii st < 10 is about 23 
percent, virtually all rays have /ii st > 2 and 93 percent of 
the rays have ni st > 5. In most of the cases where there is 
no sufficient lens plane to cause > 10 alone, there is still 
an 'almost sufficient' plane that gives rise to a magnification 
significantly larger than unity. In only 3 percent of cases is 
/i2nd > 2, and in 90 percent of all rays we find fi2nd < 1-4. 
Thus to a good approximation strong lensing can be thought 
of as being caused by individual o bjects. These result s agree 
qualitatively with the findings of IWambsganss et al.l (|2005l ) 
for the distribution of the surface mass density. 

There are, nevertheless, a few strongly lensed rays 
whose properties are due to more than one object or lens 
plane. As noted above, two or more objects at similar red- 
shift contribute significantly for a few percent of all rays. 
Fig. [9] shows that for a further few percent two or more un- 
corrected objects at different redshifts make a significant 




Figure 6. The cross section dr^/dz^d log M L for rays of type 
II V III and sources at redshift z s = 5.7 as a function of the red- 
shift z L and halo mass M L of the lens. Darker areas correspond 
to higher cross sections (on a logarithmic scale). The dashed line 
marks the mass of the largest DM halo in the Millennium Sim- 
ulation at each redshift. The solid line joins masses for which a 
spherical NFW halo of typical concentration would have Einstein 
radius 5h~ 1 kpc, approximately the resolution limit on the lens 
planes. The cross section for strong lensing in our ray sample is 
almost entirely contained between these two limits. 




M L [Mo/h] 

Figure 7. The cross section <9t : / d log M L for rays of type II V III 
to sources at 2 S = 2.1 as a function of the lensing halo mass M L . 
The solid histogram gives the direct estimate from the Millennium 
Simulation. For all other curves we have replaced each halo in the 
simulation by a spherical NFW halo with the same virial mass 
and maximal circular velocity. The heavy dotted line gives the 
result for all halos, while the heavy dashed line shows the effect of 
excluding all halos with Einstein radius < 5/i _1 kpc comoving. 
Thin lines show the result of scaling these two curves up by a 
factor of 3, so that their shape can be compared more easily to the 
direct estimate from the simulation. The cut-off in cross section 
at low halo mass in the simulation appears to correspond well to 
that induced by this simple model for the resolution limit. This 
suggests that resolution effects reduce our total cross sections by 
of order 15%. 



contribution. For the remaining rays the overall effects of 
foreground and background matter are at a much lower level. 
To demonstrate this quantitatively, we first determined the 
projected mass overdensity at the position of each ray on 
each plane. We then divided these overdensities by the crit- 
ical surface densities of the relevant planes. Finally, for each 
ray we summed the contributions from all planes to obtain 
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a 'LOS convergence' k, which - to be more precise - is the 
lensing-efnciency-weighted projection of the matter overden- 
sity along the ray. (This is equal to the convergence k in 
the single-plane and weak-lensing approximations.) We then 
performed a similar calculation for all rays of type II V III, 
with > 10 or with r > 10. In addition, when calculating 
k for these rays, we excluded the contribution of the plane 
containing the primary lens to isolate the contribution kfb 
of foreground and background matter to the lensing event. 
The cumulative distribution of the ratio of kfb to k is shown 
in Fig. 1101 For 20 percent of the rays of type II V III, addi- 
tional matter along the LOS contributes more than 10 per- 
cent to the total LOS convergence. On the other hand, for 50 
percent of the rays, there is a negative contribution kfb. Al- 
though there is a noticeable fraction of strongly lensed rays 
whose LOS convergence is enhanced by additional matter 
along the LOS, there is also a noticeable fraction of strongly 
lensed rays whose LOS convergence is decreased due to the 
lack of matter along the LOS. 

In Fig. 1111 distributions of the LOS convergence for 
z s — 5.7 are compared for all rays and for strongly lensed 
rays with the primary lens contribution removed. Although 
the distributions are very similar, small shifts are visible. 
In Fig. 1121 we show the means of these distributions as a 
function of source redshift z s . By definition, the mean LOS 
convergence of all LOS should be zero. The measured mean 
for our whole ray sample is not exactly zero because of sam- 
pling variance □ but it is much smaller than the mean for 
strongly lensed rays with the primary lens excluded. This 
demonstrates a small but measurable bias towards selecting 
directions in which matter in front or behind the primary 
lens enhances the lensing. The effect increases with increas- 
ing source redshift, is weakest for the sample of rays with 
r > 10, and is strongest for > 10. For z s — 5.7 and 

> 10, where the bias is strongest, we find an average 
contribution of 0.04 to the LOS convergence. Hence, in all 
cases the bias is small in comparison with the effect of the 
primary lens for which k ~ 1. 



4 SUMMARY 

The aim of this work has been to study the statistical distri- 
bution of the distortion of images of distant sources due to 
gravitational lensing. In particular, we have concentrated on 
estimating the cross section for rare strong-lensing events. 
Our results were obtained by shooting random rays through 
a series of lens planes creat ed from the Millennium Sim- 
ulation l|Springel et al.l 120051 ) . This is the largest iV-body 
simulation of cosmological structure formation available to- 
day. We have devised improved algorithms to make the lens- 
planes and to calculate bending angles and shear matrices on 
these planes in order to take full advantage of the very large 
volume and the high spatial and mass resolution offered by 
the simulation. 



3 One might naively expect the sampling variance to be negligi- 
ble for the 640 million rays we shot. However, the rays were not 
shot independently, but are confined to forty different patches 
of 40 X 40Mpc 2 //i 2 on each lens plane. The matter content of 
each patch is still subject to significant cosmic variance, and so, 
therefore, is the combined sampling area. 




Figure 8. The fraction of rays of type II V III (solid line), of rays 
with \fi\ > 10 (dashed line), and of rays with r > 10 (dotted line) 
to source redshift z s for which no single lens plane is able to gen- 
erate the relevant property on its own. The relative importance 
of foreground and background material clearly increases with in- 
creasing 2 s . 




Figure 9. The cumulative distribution function cdf 1 of the 
largest fii B t and of the second largest fi2nd single-plane magnifi- 
cation for all rays with a magnification \/x\ > 10 to source redshift 
2 s = 5.7. For clarity we plot 1 — cdf I (^t2nd) rather than cdf : (/i2nd) 
so that the high tail of the distribution can be compared with the 
low tail of cdf^/iist). In almost all cases the effect due to the 
primary lens is strongly dominant. 



In Sec. 13.11 we presented results for the statistical dis- 
tribution of the magnification of point sources. The dis- 
tribution is skewed with a peak at magnifications below 
unity and a tail toward high magnification. With increas- 
ing source redshift, the peak broadens and moves to lower 
magnifications, while the tail gains more weight. The mag- 
nification distribution affects the observed luminosity dis- 
tribution of astronomical standard candles. For type la su- 
pernovae, perhaps the most interesting case, magnification 
effects on the luminosity distribution are still small for cur- 
rent samples compared to the intrinsic luminosity scatter, to 
extinction and to other effects. In future high-redshift, high- 
precision surveys, however, such magnification effects may 
cause significant systematic errors, so it will be necessary 
to detect and to correct for t hem l|Dodelson fc Vallinottd 
l2006l ; lMunshi fc Vala geas 2006). In the most optimistic case, 
detailed comparison with predictions of the magnification 
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Figure 10. The cumulative distribution function cdf (kpb / /t) 
of the relative contribution kpb/k of additional matter to the 
total LOS convergence k for strongly lensed rays and for source 
redshift z s = 5.7. The dash-dotted line is for rays of type II V III, 
the dashed line for rays with \fi\ > 10, and the dotted line for rays 
with r > 10. 
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Figure 11. The probability density function pdf(SpB) of the 
LOS convergence kbf f° r strongly lensed rays and for source red- 
shift z s = 5.7, but with the primary lens contribution excluded. 
The dash-dotted line is for rays of type II V III, the dashed line for 
rays with \fi\ > 10, and the dotted line for rays with r > 10. For 
comparison, the solid line shows the corresponding distribution 
pdf(ft) for all rays irrespective of their lensing properties. A small 
but significant shift towards larger convergence is visible in the 
direction of strong-lensing events. 



distribution may help to discriminate between cosmological 
models. 

Various optical depths connected to strong lensing were 
presented in Sec. 13. 21 In particular, we estimated the fraction 
of images of sufficiently small sources that are highly magni- 
fied, have a large length-to- width ratio, or belong to multiply 
imaged sources. All the optical depths we analyse increase 
strongly with incre asing red s hift. In co mparison with ear- 
lier re sults by, e.g., iLi et al.l ([2005) and IWambsganss et al.l 
|2004h . we find a stronger evolution with source redshift, 
leading to higher optical depths for source redshifts z s > 1. 
We discussed possible reasons for this difference. 

The results we presented in Sec. 13.31 show that signifi- 
cant contributions to the strong-lensing optical depths come 
from dark matter halos with masses between 10 13 fa -1 Mq 
and 10 15 /i -1 Mq. The upper mass limit is due to the very 
rapidly decreasing abundance of more massive structures. 
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Figure 12. The mean /t mBan of the distributions shown in Fig. 1111 
as a function of source redshift z s . The line types correspond to 
those of Fig. 1111 The mean LOS convergence for all rays irrespec- 
tive of their lensing properties is non-zero only because of residual 
sampling variance. The mean convergence from foreground and 
background matter along rays containing a strong lens is signif- 
icantly non-zero, however, although still much smaller than the 
typical convergence due to the primary lens. 



This exponential decrease occurs at lower mass at higher 
redshift, and in conjunction with the lens geometry it ex- 
plains why almost all lenses are at redshifts z L < 2.5, even 
for sources with z s > 5. The lower mass limit for strong 
lensing is due primarily to the small cross sections of indi- 
vidual low-mass halos, although the spatial and mass reso- 
lution limits of the simulation itself and of our lens planes 
also play some role. Estimates based on analytic results for 
spherical NFW halos fit to the Millennium data suggest that 
these resolution effects are subdominant, and probably only 
reduce our total cross sections by of order 15%. A more im- 
portant effect on the relevant scales is our neglect of the 
baryonic mass of the central galaxies. We will come back to 
this in later work. 

We find that the mass range over which halos can cause 
strong lensing exte nds to lower masse s than those given by 
ILi et al.1 (|2005h and lDalal et~ai1 (120041 ). This difference may 
in part reflect the lower resolution of the simulations used 
by these authors, and in part the fact that they considered 
extended sources with diameters ~ 0.1-1 arcsec, while we 
assumed sufficiently small sources when calculating our cross 
sections. Halos near our lower mass limit have Einstein radii 
of order 1 arcsec and so are inefficient in producing strongly 
distorted images of sources of comparable angular extent. 

Since our set of lens planes represents the whole matter 
distribution between source and observer, we are able to 
quantify the influence of foreground and background matter 
on the frequency and the properties of strong lensing events. 
We find that such effects are quite modest. On average, the 
contribution of foreground and background material is only 
a few percent. Although we do find a bias towards excess 
foreground and background matter on strong-lensing lines- 
of-sight, the effect i s sign ificantly smaller than suggested by 
IWambsganss et all (120051 ). 

The most obvious extension to the work we have pre- 
sented here would be ray-tracing studies of the effects of lens- 
ing for realistic distributions of source properties and across 
finite size fields. This will, for example, allow direct compari- 
son with observations of massive galaxy clusters where many 
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sets of multiple images are now detected in the best cases 
l|Broadhurst et al.ll2005l ; lHalkola et alj|2006h . When galaxy 
properties from galaxy fo rmation modelling within the Mil- 



lennium Simulation (e.g . Springel et al.ll2005l; ICroton et all 



l2006l ; lBower et alJl20od ; lDe Lucia fc Blaizotll2007l ) 



are com- 
bined with such ray-tracing analyses, it will be possible to 
study whether the dark halo masses of individual cluster 
galaxies are consistent with observation, providing an addi- 
tional observational test of the hierarchical build-up of struc- 
ture predicted by the standard ACDM model. This combi- 
nation of semi-analytic simulation of galaxy formation with 
ray-tracing measures of lensing will also allow an estimate 
of how our strong-lensing cross sections should be modified 
to account for the galaxies, as well as detailed studies of pre- 
dictions for galaxy-galaxy lensing. We intend to come back 
to all these topics in future work. 
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APPENDIX A: INTEGRAL RELATIONS FOR 
THE MAGNIFICATION 

For a well behaved non-singular lens system, the numbers 
of images of type I , II and III of a giv en source satisfy 
ni — nn + Tim = 1 (Schneide r et alj|l992h . Here we briefly 
discuss an 'integral version' of this theorem for the particular 
geometry used in our work: the Multiple-Plane Approxima- 
tion with lens planes carrying a smooth and non-singular 
matter distribution that is periodic with a rectangular unit 
cell of dimensions L\ x L2. In this case, the image plane 
and source plane can both be represented by a rectangle 
P = [0, Li] x [0, £2] with periodic boundary conditions. Fur- 
thermore, the lens mapping 



9 1— > (3(G) = 9 + a(9) 



from image position 9 to source position /3 for a given 
source redshift is then smooth and non-singular with a 
smooth, non-singular, and periodic deflection angle a(9). 
Under these conditions, the inverse signed magnification 



/'• 



det(d/3/d8) satisfies the following relation: 



(Al) 



where ||P|| = L\L 2 denotes the area of the rectangle P. The 
following derivation employs integration by parts and ex- 



ploits the smoothness and periodicity of (dp/ 80): 

I d 2 9 n~ X {0) = / d 2 det 

Jv Jr 



dp 

dO 



Ae \L d6l oo ld o 2 

d0i / d0 2 

\Jo °°1 



dpi d/h 
dOi 80 

Oft. dp 2 



d<9 2 



d0i 



d0 2 



8^2 

l 86 2 

802 

90i 



=0 Jo 



9 2 =0 JO 



802 
89 2 



1 Ll 



e, =0 J 



d0i 



d0 2 0i 



d0i 



8 2 P 2 


dOu 




8 2 p 2 


80i802 


dpi 


l 2 


80! _ 


8,-- 



de 2 [p l (L 1 ,e2)-Pi(o,e 2 )} 



Li 



d9i {f3 1 (6 1 ,L 2 )-Pi(e 1 ,a)] 



89 2 

8P2(0i,L 2 ) 







802 



i2 8p 2 (L 1 ,02) 
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= L1-L2 . 

For our ray sampling method, it follows directly from 
relation (|A1|) that a representative sample of rays with ran- 
dom positions 9i (i = 1,... , N) in the image plane should 
satisfy: 

1 N 1 r 

^- 1 (80* m j/e»-\9) = l. (A2) 
For our ray sample, we find that 



- 1 < 0.003 



for all source redshifts. 

For the magnification distribution, Eq. (|A1|) implies 

that 

P df 0*) = -rrj- 



du> j fd 29 i^r 1 
j p d 2 9 K^r 1 - u(9)] 

_ 1 f F d 2 9 6[ f l' -n(9)] 



= IM 



— T \u I pdf (u ) 



= (i-2t£) iMTVdf 1 ^') • 

In practice, the optical depth 7n for images with negative 
magnification is very small. Hence, 



pdfV) « iM'r'pdf 1 ^'). 



(A3) 



Employing the fact that both probability distributions are 
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normalised, we finally find: 

l/zp 1 pdf'O) d/x « 1 , and / \fj,\ pdf S (/x) d/x « 1. 

(A4) 



APPENDIX B: EMPTY-BEAM 
MAGNIFICATION IN A FLAT UNIVERSE 

The Jacobian of the lens mapping for light propagation 
through an inho mogeneous universe is given by (see, e.g., 
ISchneiderir2006al 'l: 

Ay(0, w) = Sij 

2 r , (w-w')w' d 2 $(x(0,w'),w') , 
dw „„ „„ A k j{8,w). 



(Bl) 



Here, $ denotes the three-dimensional gravitational poten- 
tial, w the comoving line-of-sight distance and 6 the direc- 
tion of the light ray, and x the comoving transverse separa- 
tion. For an empty beam, we have A = diag(/x ml 1 1 { 2 , p m [/ 2 ). 
Furthermore, the Poisson equation in this case reads 

V 2 <f>(x(0,w),w) =-■ 3H ° 2 " m 



2a(w) ' 

where a(w) is the scale factor. 

Considering now the 1, 1-component of Eq. (IB1[) . we 

find 



m"L /2 H = i + 



dw 



2c 2 J a(w')w 



, (W - W )W -1/2, I \ 



Differentiating this expression twice, we finally obtain the 
differential equation 

d 2 / _ 3Hin m f(w) 



dw 2 2c 2 a(w) ' 

where f(w) = wy,~^^{w). This ordinary differential equa- 
tion can be easily solved numerically. 
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